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Abstract 


Th e expanded equations for torque and force on a cylindrical permanent magnet 
core in a large-gap magnetic suspension system are presented. The core is assumed to 
be uniformly magnetized, and equations are developed for two orientations of the 
magnetization vector. One orientation is parallel to the axis of symmetry, and the 
other is perpendicular to this axis. Fields and gradients produced by suspension sys- 
tem electromagnets are assumed to be calculated at a point in inertial space which 
coincides with the origin of the core axis system in its initial alignment. Fields at a 
given point in the core are defined by expanding the fields produced at the origin as a 
Taylor series. The assumption is made that the fields can be adequately defined by 
expansion up to second-order terms. Examination of the expanded equations for the 
case where the magnetization vector is perpendicular to the axis of symmetry reveals 
that some of the second-order gradient terms provide a method of generating torque 
about the axis of magnetization and therefore provide the ability to produce six- 
degree-of-freedom control. 


Introduction 

This paper develops the expanded equations for 
torque and force on a cylindrical permanent magnet core 
in a large-gap magnetic suspension system. The core is 
assumed to be uniformly magnetized, and equations are 
developed for two orientations of the magnetization vec- 
tor. One orientation is parallel to the axis of symmetry, 
and the other is perpendicular to this axis. Fields and gra- 
dients produced by suspension system electromagnets 
are assumed to be calculated at a point in inertial space 
which coincides with the origin of the core axis system in 
its initial alignment with a reference inertial axis system. 
Fields at a given point in the core are defined by expand- 
ing the fields produced at the origin as a Taylor series. 
The assumption is made that the fields can be adequately 
described by expansion up to second-order terms. The 
expansion of the fields and gradients is presented in 
appendix A. 

The equations for torques and forces on a magnetic 
core that are produced by a large-gap magnetic suspen- 
sion system have been presented and discussed in a num- 
ber of papers. For example, see references 1 through 6. 
The torques on the core are usually approximated as a 
function of the external or applied fields at the centroid 
of the core, and the forces on the core are usually approx- 
imated as a function of the gradients of the applied fields 
at the centroid. It is generally assumed that terms that are 
a function of second-order or higher gradients of the 
applied fields at the centroid can be neglected. In practi- 
cal applications that involve large-gap magnetic suspen- 
sion systems, these assumptions have proven to be valid. 
For an axisymmetric core, such as a cylinder, it can be 
shown that if the direction of magnetization is along the 
axis of symmetry, then the torque about that axis, pro- 
duced by the applied fields and gradients of the applied 
fields, is always zero (ref. 3). Various methods of over- 


coming this constraint, which include shaping the core 
and using nonuniform three-dimensional magnetization, 
are discussed in references 4 and 7. However, examina- 
tion of the expanded equations for a cylindrical core 
reveals that for the case of uniform magnetization per- 
pendicular to the axis of symmetry, some of the second- 
order gradient terms provide a method of generating 
torque about the axis of magnetization and therefore the 
ability to produce six-degree-of-freedom control (ref. 8). 
For completeness, all gradient terms for expansion of the 
fields up to second order are presented. 

Finally, instead of developing the torque and force 
equations from a set of governing equations that are a 
function of core volume, core magnetization vector, and 
suspension system fields and gradients, appendix B pre- 
sents a development that begins at a more fundamental 
level in an attempt to provide better insight into the ori- 
gin of these equations than is commonly available in the 
literature. 

Symbols 

A area, m 2 

a radius of core, m 

B magnetic flux density vector, T 

B expanded magnetic flux density vector, T 

[dB] matrix of field gradients, T/m 

[dB ] matrix of expanded field gradients, T/m 

F total force vector on core, N 

5F force vector on incremental volume of core, N 

I coil current vector, A 

l length of core, m 

m magnetic moment vector, A-m 2 


8m magnetic moment vector of dipole with incre- 

mental volume, A-m 2 

M magnetization vector. A/m 

Q m pole strength, A-m 

r position vector, m 

T total torque vector on core, N-m 

8T torque vector on incremental volume of core, 

N-m 

ST torque on incremental volume of core about 

core origin, N-m 

T m inertial coordinate to suspended-element coor- 

dinate vector- transformation matrix 

U potential energy 

v core volume, m 

W work 

5v incremental volume, m 

x, y, z coordinates in orthogonal axis system, m 
0 Euler orientation for 3, 2, 1 rotation sequence, 

rad 

V gradient operator 


Subscripts: 

x, y, z components along x -, y- t z- axes, respectively 
ij partial derivative of i component in ./‘-direction 

(ij)k partial derivative of ij partial derivative in 

^-direction 

Matrix notations: 

[ ] r transpose of matrix 

LJ row vector 

A bar over a symbol indicates that it is referenced to 
suspended-element coordinates. 

Magnetic Torques and Forces 

The torques and forces on a cylindrical permanent 
magnet core are developed in this section by integrating 
the equations for torques and forces on an incremental 
volume of the core with magnetic moment M8v over the 
core volume. These equations are developed in appendix 
B. Figure 1 shows the cylindrical core and the core coor- 
dinate system. The core coordinate system consists of a 
set of orthogonal x , z body-fixed axes that are initially 

aligned with a set of orthogonal x-, y-, z-axes fixed in 
inertial space. In order to define the fields and gradients 
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at any point in the core, the fields and gradients at the ori- 
gin of the core axis system are expanded as a Taylor 
series. It is assumed that the fields can be adequately 
described by expansion up to second-order terms. The 
expanded fields and gradients in both inertial and core 
coordinates are presented in appendix A. For simplicity 
in developing the equations in this section, relative 
motion between the core and the reference inertial coor- 
dinate system is assumed to be zero. This assumption 
removes the requirement to transform between the iner- 
tial and core coordinate systems and eliminates a signifi- 
cant number of components which are small relative to 
the fundamental terms in the equations when small-angle 
assumptions are used. In particular, the transformation of 
second-order gradient terms from inertial to core coordi- 
nates is very complicated, as illustrated by equa- 
tion (A 14). The torque on an incremental volume of the 
core, about the core origin, can be written as 

ST = 8T + (rx8F) (1) 

where 8T and 8F are the torque and force on the incre- 
mental volume due to the field at that location, and 


r - 


is the position vector of the incremental volume 
(fig. 2). The total torque on the core can be written as 


T = J ST = J [8T + (r x 8F)] (2) 

V V 

where the integration is over the core. Substituting equa- 
tions (B20) and (B21) results in 


T - J {(M x B) + [r x (M • V)B] } dv 

V 

The total force on the core can be written as 
F = J (M* V)B dv 

V 

The term (M • V)B can be written as (ref. 5) 
(M • V)B = [dB]M 


( 3 ) 


( 4 ) 


( 5 ) 
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where 


_3b] “ 


Bjtx 

By* ®yy ®y; 
B zjc ®zy ®z i 


( 6 ) 


and 


M = 


M .V 

M- 


(7) 


In equation (6) the notation fy = df/dj is used. 

Magnetization Along Axis of Symmetry 

For orientation of the magnetization vector along the 
axis of symmetry (x-axis) of thepermanent magnet core, 
the only nonzero term in M is M-. Expanding 
equation (5) results in 


[3B]M = M- 




( 8 ) 


Substituting equation (8) into the second part of equation 
(3) and expanding results in 


r x (M • V)B = M- 


(-Bjyj + B rz y) 
(B^Z — B,. z x) 
[~B xx y + B tyx) 


(9) 


Expanding the first part of equation (3) results in 


MxB = M ; 


0 

-B* 

B, 


( 10 ) 


The components of equation (3) become 

T^ = 0 + M- J (B„y-B xy z) 


dv 


(id 


T- = -M- J B z dv + M i J (B xx z - B xz x) dv (12) 
Tj = J By dv + J (BjtyX - B xx y) dv (13) 

V V 

Evaluating J B xz y dv first 

V 

J B, z y dv = B xz \ydv + B {xx)z J xy dv 

V V v 

+ B (x»z\y 2 dv+B (szu\* ydv (14) 


where B xz has been expanded by using a Taylor series, 
as detailed in appendix A, and the notation 
f = d(df t /dj)/dk has been used. All the integrals 
involving first-order terms in equation (14) are zero. 
Evaluating J y 2 dv yields 


f ? dv = J t. dldi 


- 1/2 


( a /4)nl 


(15) 


where y,(z) = y 2 (z) = Ja*-z 2 , a is the radius of the 
permanent magnet core, and l is the length (fig. 1). Since 
the area of the face of the permanent magnet core is 
A = n o' and the volume is v = Al, equation (15) 
reduces to 

(16) 


(a 4 /4)tt/ = (a 2 /4)v 


(17) 


m 


Substituting in equation (14) results in 

J B^y = {a / 4)vB (jty)z 

V’ 

Evaluating | B xy z dv in a similar manner results 

V 

J B xv .z dv = (a 2 / 4)vB (x> , )z (18) 


which is equal to equation (17). Therefore the torques 
about the x-axis due to second-order gradients cancel 
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out and T- = 0 as expected. Going next to equa- 
tion (12), 


Finally, noting that B (ij)k = B (ik)j = B (Jk)i ... and collect- 
ing terms, the components of torque become 


V 

+ (‘’ 2 / ( 19 ) 

and 

\(B xx i-B xz x)dv = v[{a/A)-(l 2 /\2)}B (xx)z (20) 

V 

Substituting into equation (12) results in 

T- = - V M,[B z + U 2 /24)B (zx)x + ( a 2 /S)B (zy)y 

+ (« 2 /8)B (zz)z + vM-|~(a 2 /4) 

-( / 2 /12)lB (21) 

J ( xx)z 

Continuing to equation (13), 


T i = 0 (25) 

T- = -vM-B z - vM-[(« 2 /4)-(/ 2 /8)]B (jrjc)j 

- vM-(« 2 /8)(B (y> , )z + B (zz) J (26) 

T -z = ‘'M-B^+vM- t( /2/8 )-(« 2 /4)]B (jr ^ 

- vM i( « 2 /8)(B (w)y + B (vz)z ) (27) 

The force on the core, equation (4), can be evaluated 
in a similar manner. The components of equation (4) are 

F i = M jB*^ (28) 

V 

Fy = M_jB jy Jw (29) 

V 

F- = M_J§ x ,</v (30) 


jB y dv= v[B y + (l 2 /24)B {yx)x + (a 2 /8)B (yy)y 


Expanding the integral of equation (28) results in 


V 

+ (« 2 /8)B (yz) ^ 

(22) 

J® xxdv = J( B xx + B («)^ + B (jtJt)z z 

V V 




+ B (**)x* dv 

(31) 

and 


Since the integrals containing first-order terms in 

x, y. 



and z are zero, equation (31) reduces to 


J(B iy *-B„y)rfv = v[(/ / 12) - (a / 4) ]B (jtx jy 

(23) 




F- = M_ f B dv = vM B 

x X J XX x XX 

(32) 

Substituting into equation (13) results in 


V 




Evaluating equations (29) and (30) results in 


Tj = vM s [B > + (/ 2 /24)B (y , )y + (fl 2 /8)B (yy)y 


F » - vM.B v 

(33) 

+ (» 2 /8)B (b) J + .M i [(/ 2 /12) 


and 



(24) 

F, = vM_B 

Z X xz 

(34) 
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Magnetization Perpendicular to Axis of The components of the force (eq. (4)) using equa- 

Symmetry tion (35) becomes 


For orientation of the magnetization vector perpen- 
dicular to the axis of symmetry ( v- axis), the only nonzero 
term in M is Equation (5) becomes 


[3B]M = M- 


xz 


yz 


B, 


(35) 


F i = M j®« rfV = VM i B XZ 

V 

(44) 

F y — M z J By Z dv = vM x B n 

V 

(45) 

F. = M J§ zz rfv = vM,B,, 

(46) 


V 


and the second part of equation (3) becomes 


r x (M • V)B = M: 


(-B >z z + B, z y) 

(BjczZ-Bjj*) 

(-B, z y+B vz x) 


The first part of equation (3) becomes 


MxB = M- 


-B, 


B, 

0 


(36) 


(37) 


Substituting (36) and (37) into (3) results in 
T- = -MjB y dv + Mj(B zz y - B yl z)dv 


(38) 


T- = M- J B x dv + M- J (B 


«*-B zz x)dv (39) 


T- = 0 + M- J (By Z x - B xz y)dv 


(40) 


Evaluating the integrals as before and collecting terms 
results in 


T* = -vM-B y -vM-(/ 2 /24)B (jtx)y . 

_vM,(a 2 /8)(B ( ^ )> , + B (> , z)z ) (41) 

T- = vM-B x +vM-[(3a 2 /8) -(/ 2 /12 )]B (xz)z 

+ vM-(/ 2 /24)B (xx)x + vM> 2 /8)B ( ^ )y (42) 

T- = vM.[(/ 2 / 12) - (a 2 /4) ]B (xy)z (43) 


Discussion of Results 

Examination of equations (25) to (27), (32) to (34), 
(41) to (43), and (44) to (46) reveals that, for expansion 
of the applied fields up to second-order terms, no cou- 
pling exists between force and torque components. As 
stated earlier, it is generally assumed that the higher 
order torque terms, which are functions of second order 
gradients, can be neglected. However, equation (43) indi- 
cates that for magnetization perpendicular to the axis of 
symmetry, torque about the axis of magnetization can be 
generated by controlling a higher order term directly, 
thus allowing the core to be controlled in six degrees of 
freedom. For a cylindrical permanent magnet core mag- 
netized along the axis of symmetry, from equation (25), 
only five-degree-of-freedom control is possible. 

Concluding Remarks 

This paper has developed the expanded equations for 
torque and force on a cylindrical permanent magnet core 
in a large-gap magnetic suspension system. The core was 
assumed to be uniformly magnetized, and equations were 
developed for two orientations of the magnetization 
vector. One orientation was parallel to the axis of sym- 
metry of the core and the other was perpendicular to this 
axis. It is generally assumed that terms that are a function 
of second-order or higher gradients of the applied fields 
can be neglected. In practical applications involving 
large-gap magnetic suspension systems, these assump- 
tions have proven to be valid. However, in the case 
where the magnetization vector is perpendicular to the 
axis of symmetry of the core, the expanded equations 
indicate that torque about the magnetization vector can 
be produced by controlling a second-order gradient 
directly. This case allows the core to be controlled in six 
degrees of freedom whereas a cylindrical permanent 
magnet core magnetized along its axis of symmetry can 
be controlled only in five degrees of freedom. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
October 24, 1996 
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Appendix A 

Expansion of Fields and Gradients About the 
Nominal Operating Point of a Cylindrical Per- 
manent Magnet Core 

In appendix A the fields and gradients produced by 
the suspension system electromagnets are expanded by 
using a Taylor series about the initial suspension point of 
the permanent magnet core. The assumption is made that 
the fields can be adequately described by expansion up to 
second-order terms. Figure 1 shows the cylindrical core 
and core coordinate system. The core coordinate system 
consists of a set of orthogonal x , y, z body-fixed axes 
that define the motion of the core with respect to an 
orthogonal x , y, z system fixed in inertial space. The core 
coordinate system is initially aligned with the x, y, z sys- 
tem. The transformation from inertial coordinates to core 
coordinates is given by 


X 

= IT J y (Al) 

A 

where [ T m i is the orthogonal transformation matrix for 
a 3, 2, 1 (z, y, x) Euler rotation sequence and is defined as 



x 

where r = y and V is the gradient operator. Using 

z 

compact notation, each element of B in equation (A3) 
can be written as 

dB, T 3 2 B,. 

Bj < = B . + a7 r + (1/2)r ~T r (A4) 

dr 

where 

dB. dB. dB, dB. 

dF = Lai a7 aFJ (A5) 

and 


d(GB,./dx) 

d(dB /dx) 

d(dB,/dx) 

dx 

dy 

dz 

3(dB /dy) 

d(dB/dy) 

d(dB -/dy) 

dx 

dy 

dz 

d(dB /dz) 

G(GB-/dz) 

d(dB -/dz) 

dx 

dy 

dz 


Using the notation f (j = df ( /dj and f k = 
d(df i /dj)/dk, equations (A5) and (A6) can be written 
as 


ce,c0 v s0,c0 y -s0 v 

f (cG^G^G^-sG^GJ (sG^BySG^ + cG.cG^) cO y sQ x 
(cG^sG^cG^ + sB^sG^) (sB^sG^.cG^*- e0,s0J cG^cB^ 

(A2) 

where sin has been shortened to s, cos has been shortened 
to c, and 0^, 0 y , and 0^ are angles of rotation about the 
Z-, y-, and x-axes, respectively. The field B and gradients 
of B produced by the suspension system electromagnets, 
which are fixed in the inertial frame, are calculated at the 
origin of the x, y, z system. 

Expanding B about the origin of the x, y, z system as a 
Taylor series, up to second order, results in 



B = B + (r.V)B + (l/2)(r«V) 2 B (A3) 


9(3B/d/) 

=r = B,., B, ... B,... (A10) 

dr (‘J) x ( <j)y (y)z v ' 
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The expanded fields can be expressed in core Transforming back into core coordinates, 
coordinates as 


B = B + (r • V)B + (l/2)(r • V) B (All) 


where r = y is the displacement in core coordinates, 

z_ 

B = [T m ]B, and V = [T m ]V. Since, from equa- 
tion (Al), 


T- 

r = [TJ r 


(A 12) 


each element of B can be expanded in inertial coordi- 
nates by substituting equation (A 12) into equation (A4), 


B, = B | . + SrTJ r f + (l/2)f T [T w li4 , [T Jll l r f (A13) 


j ^ B t *. 


B, + ¥ Ul ^ (1/2)r fTJ— ^TJ i 


H rT 1 vB,, T _ T v B T — 

B = [T « 1 B v + -^[TJ r f + (l/2)r [TJ-^fTJ r 


2 

0B„ t -T 0 B, r 

B * + 3F rT "i r r + 0/2)F [TJ-f[T m l r 


(A 14) 

The expansion of equation (A 14) can be simplified by 
using small-angle assumptions (ref. 6). Under small- 
angle assumptions, cos0 = 1, sin0 = 0, and products 
of angles are neglected. The transformation matrix [T m ] 
then becomes 


6 Z - 6 , 


L e >- - e * 1 
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Appendix B 

Torques and Forces on a Magnetic Dipole 
With Incremental Volume 

The torques and forces on a magnetic dipole in a 
steady magnetic field are identical to those on an infini- 
tesimal current loop with the same magnetic moment 
(ref. 9). Therefore, the equations for torque and force on 
an infinitesimal current loop will be developed first by 
using the fundamental relationship for the force on a 
current-carrying-conductor element in a uniform, steady 
magnetic field. For a discussion of magnetic dipoles and 
infinitesimal current loops, see references 9 and 10. 

Infinitesimal Current Loop 

Consider a plane loop of conductor with steady cur- 
rent I located in the external, uniform, steady magnetic 
field B (fig. Bl). In this region VxB = V*B = 0. 
The force on an element d\ of the conductor is given by 
the fundamental relationship (obtained from the Lorentz 
force law) 

d¥ = I dl x B (Bl) 

where dF is a vector indicating magnitude and direction 
of force on the conductor element; I is the scalar magni- 
tude of the current in the conductor element; d\ is a vec- 
tor whose magnitude equals the length of the conductor 
element and whose direction is in the positive direction 
of the current; and B is a vector indicating magnitude and 
direction of the flux density of the external field compo- 
nent. The torque on the loop can be written as 

T = I<j>[r x(<flxB)] (B2) 

where r is the position vector of d\ and the integration is 
around the loop. By using the identity 

rx(JlxB) = d\{v • B) - B(r • d\) (B3) 

equation (B2) can be written as 

T = I[£(r • B)dl - B^r • dl] (B4) 

Using Stokes’s theorem and a related result (ref. 9, 
p. 289), the line integrals in equation (B4) can be trans- 
formed into surface integrals resulting in 

T = ij J [d A x V(r • B)] -B |(V x r) • dA 1 (B5) 

^ S S ' 


where dA is a vector whose magnitude is a differential 
area and whose direction is normal to the plane of the 
current loop in the sense of the right-hand rule relative to 
the direction of current flow, V is the gradient operator, 
and the integrals are over the surface that is defined 
by the conductor loop. Since V x r is zero and 
V(r • B) = B for constant B, equation (B5) simplifies 
to 

T = I J (dA x B) (B6) 

S 

Taking the integral results in 

T = IA x B (B7) 

An infinitesimal current loop can be defined by letting A 
go toward zero and I go toward infinity, keeping the 
product IA finite. For an infinitesimal current loop, the 
requirement that B be uniform no longer exists. The 
product IA is called the magnetic moment of the loop 
and is designated by the letter m. Therefore equation 
(B7) becomes 

T = mxB (B8) 

The torque T acts on the infinitesimal current loop in a 
direction to align the magnetic moment m with the exter- 
nal field B. If m and B are misaligned by the angle 0, the 
magnitude of the torque is 

T = mBsinG (B9) 

To increase 0 by the amount J0, work dW must be done 
against the torque T resulting in an increase in potential 
energy dU: 

dU - dW = T dQ = mBsinGdG (BIO) 

The potential energy of an infinitesimal current loop in 
an external magnetic field can then be obtained by inte- 
grating equation (BIO): 

U = -mBcosG = -m • B (B 11) 

where the constant of integration is chosen to be zero 
when m is perpendicular to B. The force on the infinites- 
imal current loop can be obtained from equation (B 1 1 ). If 
an external force F displaces the infinitesimal current 
loop by the infinitesimal distance dr, then the work done 
dW will be equal to a decrease in potential energy, -dU\ 

dW = F^dr = -dU = -VU •dr (B12) 

Therefore 

F = -VU = V(m • B) (B 13) 
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The right-hand side of equation (B13) can be expanded 
as 

V(m • B) = mx(V x B) + (m • V)B 

+ Bx(Vxm) + (B • V)m (B14) 

Since (V • B), (V x B), and (V x m) are zero, equa- 
tion (B13) can be written in the form 

F = (m • V)B (B15) 

This form is generally used in the development of the 
equations for large-gap magnetic suspension systems 
(refs. 1 through 6). 

Magnetic Dipole With Incremental Volume 

The magnetic moment of a permanent magnet dipole 
with north and south poles separated by length / and with 
pole strength Q m is defined as 

m = Q m l (B16) 

The magnetic moment m is a vector pointing from the 
south pole to the north pole. In the case of an actual mag- 
net, Q m and / may be indefinite but m can be determined 
and is sufficient to specify the fields of the magnet at a 
large distance from it. At large distances, a magnetic 
dipole with magnetic moment Q m l can be treated the 
same as an infinitesimal current loop with magnetic 
moment I A and is identical in effect if Q m l - I A. There- 
fore, in a steady magnetic field B the equations for 
torques and forces on a magnetic dipole with magnetic 
moment m are the same as equations (B8) and (B15). 

In theory, it can be assumed that a permanent magnet 
of a given volume v consists of a large number of uni- 
formly distributed permanent magnet dipoles with incre- 
mental volumes 8 v which are oriented in the same 


direction. The magnetic moment 5m of a given dipole 
with incremental volume 8v can be conveniently 
described by a quantity called the magnetization M, 
which is defined as the magnetic moment per unit vol- 
ume. That is, 

M = 8m/5v (B17) 

The total magnetic moment m for a given permanent 
magnet can then be written as 

m = J M dv (B18) 

V 

where the integration is over the volume of the perma- 
nent magnet. Magnetization is also a vector and has the 
same direction as m. If the permanent magnet is uni- 
formly magnetized, that is, M is constant over the vol- 
ume of the permanent magnet, then 

m = Mv (B19) 

For a discussion of magnetic dipoles and magnetization, 
see reference 10. 

The torques and forces on an incremental volume of 
permanent magnet material, in terms of the magnetiza- 
tion M, can then be written as 

8T = (M x B)5v (B20) 

and 

5F = (M • V)B8v (B21) 

from equations (B8), (B15), and (B17). 
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